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Abstract. Plateaus can be observed in the zero-temperature magnetization curve of 
quantum spin systems at rational values of the magnetization. In one dimension, the 
appearance of a plateau is controlled by a quantization condition for the magnetization 
which involves the length of the local spin and the volume of a translational unit cell 
of the ground state. We discuss examples of geometrically frustrated quantum spin 
systems with large (in general unbounded) periodicities of spontaneous breaking of 
translational symmetry in the ground state. 

In two dimensions, we discuss the square, triangular and Kagome lattices using 
exact diagonalization (ED) for up to = 40 sites. For the spin-1/2 XX Z model on 
the triangular lattice we study the nature and stability region of a plateau at one third 
of the saturation magnetization. The Kagome lattice gives rise to particularly rich 
behaviour with several plateaus in the magnetization curve and a jump due to local 
magnon excitations just below saturation. 
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1. Introduction 

Geometrically frustrated quantum spin systems constitute a class of systems exhibiting 
interesting quantum phenomena like different (unusual) quantum phases (see e.g. [1,2] 
for reviews of two-dimensional models). An external magnetic field may further enhance 
this frustration since it competes with the antiparallel alignment of spins favoured 
by antiferromagnetic exchange. This leads to interesting quantum phenomena in the 
high-field magnetization process of geometrically frustrated quantum spin systems. In 
particular, a spin gap may be opened by the external magnetic field, giving rise to a 
plateau in the T = magnetization curve. Although no complete review of the subject 
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exists so far, different aspects are summarized e.g. in [1-5]. In this article we will 
illustrate the status of the field by pointing out and presenting selected recent and new 
results. 

Specifically we consider the XX Z model in an external magnetic field h 

where the Si are spin-S* operators at site i and the Jij are the exchange constants 
between pairs of sites {i,j) connected in the lattice topology under consideration. The 
special case of (P) where the XXZ anisotropy satisfies A = 1 corresponds to the 
Heisenberg model. Below we will concentrate on the extreme quantum limit S = 1/2 
and antiferromagnetic coupling Jij > although some of the results can be generalized 
e.g. to S > 1/2. 

An important observable is the magnetization 

which we normalize to saturation value (M) = 1 (here is the total number of spins in 
the system). Note that the magnetization (0) is a conserved quantity for the Hamiltonian 
This is technically useful for computing the magnetization curve e.g. by exact 
diagonalization (ED) using the Lanczos method. 

2. One dimension 

In one dimension, Oshikawa, Yamanaka and Affleck proposed the following condition 
on the magnetization (M) at a plateau in the T = magnetization curve [6] 

SV{l-{M))eZ. (3) 

Here 5* is the size of the local spin and V is the number of spins in a translational unit 
cell of the ground state. The condition can be understood easily starting from a 
limit where the system decouples into clusters with V sites (see e.g. [7]) and in this case 
V is the number of spins in a unit cell of the Hamiltonian. However, in particular in 
frustrated spin systems one finds that translational symmetry is spontaneously broken 
in the ground state (c.f. examples in the following sections) and then V is larger than 
(namely an integer multiple of) the unit cell of the Hamiltonian. Nevertheless, the 
quantization condition implies that the magnetization (M) always has a rational 
value on a plateau. 

2.1. Frustrated S = 1/2 Heisenberg chain 

Let us start by discussing some selected aspects of the = 1/2 frustrated Heisenberg 
chain (see inset of Fig. Q). Majumdar and Ghosh noticed more than 30 years ago [8-10] 
that one can explicitly write down the ground state if the next-nearest neighbour 
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Figure 1. Inset: Heisenberg chain with, nearest and next-nearest neighbor coupling 
Ji and J2, respectively. Points denote the locations of S* — 1/2 spins. Main panel: 
Magnetization curve for Ji — 1, J2 — 0.8. The thin dashed and solid curves were 
obtained by ED of rings (periodic boundary conditions) with L = 12 (dashed), 24 
(dotted) [12] and 36 sites (full); the bold solid curve was obtained by DMRG for 
L = 192 sites with open boundary conditions [13]. 

coupling J2 and the nearest neighbour couphng Ji are related by J2 = Ji/2. This 
so-called 'Majumdar-Ghosh' state is a two- fold degenerate dimer product state and has 
become a paradigm for the frustration-induced opening of a spin gap accompanied by 
spontaneous breaking of translational symmetry by a period two. 

In the presence of a magnetic field, the spin gap leads to an (M) = plateau in the 
magnetization curve. Starting with [11], many authors have studied the magnetization 
process of this model by ED (see e.g. [12] and references therein) and no plateaus with 
(M) 7^ were found. Finally, a recent density-matrix renormalization group (DMRG) 
study of the 5 = 1/2 frustrated Heisenberg chain exhibited a clear plateau at (M) = 1/3 
in the region of strong frustration Ji ^ J2 [13], as is illustrated by the main panel 
of Fig. ^ which reproduces the magnetization curve for J2 = 0.8 Ji (the step on the 
(M) = 1/3 plateau in the DMRG curve is due to the open boundaries). The state of 
this (M) = 1/3 plateau spontaneously breaks translational symmetry by a period three. 
Two points are worthwhile noticing here. First, this period three is higher than the 
period two appearing in the Majumdar-Ghosh state [8-10] at h = 0. Second, the very 
presence of this plateau was missed for a long time. 

With hindsight, this (M) = 1/3 plateau can be observed also in the ED results for 
L = 12 and 24 sites and periodic boundary conditions (reproduced for J2 = 0.8 Ji in 
Fig. Rafter [12] - this figure also includes a new L = 36 curve). It may be instructive to 
consider the reasons for missing it nevertheless. First, the excited states at h = [14] 
and correspondingly the ground states in a magnetic field h ^ have in general 
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Figure 2. Inset: Orthogonal dimer chain. Points denote the locations of S — 1/2 
spins. Main panel: Magnetization curve for J' — 0.7 J [19]. Note the long (probably 
infinite) sequence of plateaus for (Af) < 1/2. 

incommensurate momenta [12] which leads to oscillating finite-size effects. Second, one 
should concentrate on system sizes that can accommodate the known antiferromagnetic 
ground states in the limit of two decoupled chains Ji — > 0, implying that L should be a 
multiple of 4. In addition, only those values of the magnetization can be studied which 
are realized for a given system size. In particular, (M) = 1/3 is realized only if L is a 
multiple of 3. These considerations restrict the system sizes that should be used for the 
discussion of an (M) = 1/3 plateau in the frustrated Heisenberg chain to multiples of 
12. The first cases are L = 12, 24 and 36. These are used in Fig. ^ and are in good 
agreement with the DMRG results of [13]. 

2.2. Orthogonal dimer chain 

Now we wish to draw the reader's attention to the 5 = 1/2 orthogonal dimer chain 
(also known as dimer-plaquette chain) [15, 16] shown in the inset of Fig. |21 a further 
one-dimensional model where some remarkable properties were missed originally. This 
model attracted renewed attention recently since it can be regarded as a one-dimensional 
version of the two-dimensional Shastry-Sutherland model [17]. 

A first study [18] of the magnetization process of the S" = 1/2 orthogonal dimer 
chain found plateaus with (M) = 1/4 and 1/2. However, this study missed further 
plateaus and thus the remarkable property of the orthogonal dimer chain that it gives 
rise to a magnetization curve with (probably infinitely) many plateaus for some values 
of the parameters [19]. At least some steps in the derivation can be carried out exactly 
exploiting the special property of the orthogonal dimer chain that the total spin on each 
vertical dimer is conserved. A plateau-state can then be characterized as follows [19]: 
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k consecutive vertical dimers form triplets {S = 1), separated by singlets {S = 0). 
Covering the chain periodically with such 'fragments' of length k each in the sector 
= k yields 

At J' = 0.7J > and for (M) < 1/2 only precisely these states consisting of fragments 
of length k with = k appear as ground states in a magnetic field [19], whereas 
for (M) > 1/2 the ground state is in a sector where all vertical dimers form triplets 
S = 1 (one difference between the present Fig. |21 and the original Fig. 3 of [19] is in 
the system sizes that have been used for extrapolating the magnetization curve in the 
region (M) > 1/2). As long as (M) < 1/2, fragments of increasing length become 
ground states with increasing field, leading to a magnetization curve with infinitely 
many plateaus at the magnetization values given by Q and jumps in between, compare 
the main panel of Fig. |2l This infinite series accumulates at (M) = 1/2, where one finds 
a further pronounced plateau before a smooth transition to saturation follows. 

The infinite sequence of plateaus arises in the orthogonal dimer chain because the 
ground states manifestly break translational symmetry with arbitrarily long periodicity 
A; + 1, implying that at least in one dimension there is no general upper bound on the 
possible periods for spontaneous breaking of translational symmetry. 



3. Two dimensions 



To the best of our knowledge there are no rigorous arguments why the condition (jH]) 
should hold in dimensions higher than one. Furthermore, starting from two dimensions 
one may have true long-range order in the spin components at T = 0. Accordingly, an 
external magnetic field can induce transitions between different ordered and disordered 
states. 

In one dimension, a plateau at a given (M) already requires a detailed analysis of 
the specific model. And, even more so the presence and nature of a plateau state must 
be discussed case by case in two dimensions. Some selected examples are presented in 
the following sections. 



3.1. Square lattice 

Let us start with a brief discussion of the magnetization process of the square lattice 
which is a non-frustrated model and well understood. Fig. El shows the magnetization 
curve of the 5 = 1/2 square-lattice Heisenberg model (A = 1) obtained by different 
approaches. First, the thin full line shows the result obtained by ED for a finite lattice 
with iV = 40 sites (see also [20,21] for earlier ED studies - in the particular case N = 40, 
the largest dimension to be treated was 554 596 326 and occurred in the sector = 3). 
The full bold line denotes an extrapolation of the ED data to the thermodynamic limit 
which is obtained by connecting the midpoints of the finite-size steps at the largest 
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Figure 3. Magnetization curve of the S ~ 1/2 Heisenberg antiferromagnet on 
the square lattice. The thin sohd Une is for = 40 sites, the fuU bold line is an 
extrapolation to the thermodynamic limit. Also shown are a second-order spinwave 
result [22] (dashed bold line) and QMC results [23] (diamonds). 

available system size. One observes a smooth magnetization curve with no peculiar 
features (in particular no plateaus) for |(M)| < 1. Note that close to saturation the 
extrapolated curve includes data at large system sizes, which is not shown explicitly 
in Fig. ini More precisely, for (M) > 0.84375, the curve is based exclusively on finite 
lattices with at least 8x8 sites. The high-field part of the magnetization curve is 
therefore particularly well controlled by ED. 

Second, a second-order spinwave result [22] is also available and shown by the bold 
dashed line in Fig. El Third, the magnetization process of the square lattice can also 
be studied by quantum Monte Carlo (QMC) since this lattice is not frustrated. The 
diamonds in Fig. El show available stochastic-series-expansion QMC results [23] . 

The quantitive differences of the results of all three approaches are small, i.e. each 
approach yields a good description of the S' = 1/2 square lattice. Since the spinwave 
approach [22] is based on a Neel state, we may therefore conclude that Neel order prevails 
in the transverse components for | (M ) | < 1 (see also [20] for a discussion from the point 
of view of ED). 

3.2. Triangular lattice 

The 5=1/2 XX Z model on the triangular lattice is among the first models whose 
magnetization process was studied by ED [24]. These early studies already found a 
plateau with (M) = 1/3, at least for Ising-like anisotropies A > 1. Due to the restriction 
to at most 21 sites, it was first not completely clear whether the plateau persists in the 
isotropic regime A ^ 1. The magnetization process of the isotropic model was analyzed 
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Figure 4. Magnetization curve of the S = 1/2 Heisenberg antiferromagnet on the 
triangular lattice. The thin dashed and solid line are for N — 27 and = 36 sites, 
respectively. The bold line is an extrapolation to the thermodynamic limit. 

further using spinwave theory [25]. This study provided evidence that the (M) = 1/3 
plateau exists also at A = 1 and estimates for its boundaries were obtained. 

Fig. m shows the magnetization curves obtained by ED for A = 1 on finite lattices 
with N = 27 and 36 sites (thin lines). Both curves exhibit a clear plateau at (M) = 1/3 
in an otherwise smooth magnetization curve. These two curves overlap with previous 
ED results [21,26]. Note, however, that for = 36 only the sector k = was studied 
previously for (M) < 1/3 [26], but it turns out that the ground state is not always in 
this sector. Inclusion of sectors with k for = 36, a different shape of the N = 27 
system and a previously incomplete = 36 curve lead to small differences of our Fig. |3] 
as compared to Fig. 4 of [21]. The first-order spinwave results for the magnetic fields at 
the lower and the upper boundaries of the (M) = 1/3 plateau [25] are smaller by about 
0.13J (lower boundary) and O.OIJ (upper boundary) than the ED results presented here 
for = 36 and 5 = 1/2. 

The full bold line in Fig. H] denotes an extrapolation of the ED data to the 
thermodynamic limit which is obtained by connecting the midpoints of the finite-size 
steps at the largest available system size (except for the boundaries of the (M) = 1/3 
plateau where corners were used) . Close to saturation this includes again bigger system 
sizes than those explicitly shown in Fig. |3 

The state of the (M) = 1/3 plateau can be easily understood in the Ising limit 
A ^ 1 [21,27]. Quantum fluctuations are completely suppressed in the limit A ^ oo 
and the (M) = 1/3 state is a classical state where all spins on two of the three sublattices 
of the triangular lattice point up and all spins on the third sublattice point down, as is 
sketched in the inset of Fig. El This state corresponds to an ordered coUinear spin 
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Figure 5. Overlap of the wave function |full XXZ) of the full 5 = 1/2 XX Z model 
on the triangular lattice at (M) = 1/3 with the corresponding state |Ising) of the 
Ising model as a function of the XXZ anisotropy A. Lines are for finite lattices with 
TV = 18 (dotted), 27 (dashed) and 36 (full) sites. The inset shows one representative 
state |Ising) for the three degenerate ground states of the Ising model at (M) = 1/3. 

configuration. It breaks translational symmetry by a period three and accordingly 
is threefold degenerate. One can then use perturbation theory in 1/A to study the 
(M) = 1/3 plateau of the XXZ model [21]. Here we numerically compute the overlap 
of the appropriate linear combination | Ising) of the three Ising states and the (M) = 1/3 
wave function |full XXZ) of the full XXZ model with S = 1/2 in order to show that 
this description remains qualitatively valid even in the isotropic region A ^ 1. Fig. El 
shows results for the overlap | (full XXZ | Ising) | on finite lattices with N = 18, 27 and 36 
sites (for X = 36 we have computed the overlaps by considering only the wave function 
within the /c-subspace that contains the ground state for A > 1). One observes that 
this overlap remains large even in the vicinity of the Heisenberg model A = 1 and drops 
sharply around A = 0.75. Accordingly, this (M) = 1/3 plateau state is a stable phase 
for A > Ac with Ac ~ 0.75. Remarkably, this estimate for Ac is close to the estimate 
Ac ~ 0.85 for the vanishing point of the (M) = 1/3 plateau which was obtained in [21] 
by comparing the plateau widths on 3 x 6 and 3x9 systems. Note that the location 
of the sharp drop in Fig. El is almost independent of system size although the absolute 
values of the overlaps |(full XXZ | Ising) | do depend on X. We can thus locate the point 
Ac where the (M) = 1/3 plateau disappears quite accurately at Ac = 0.76 ± 0.03. 

3.3. Kagome lattice 

Finally, Fig. IHl shows magnetization curves for the Kagome lattice (see inset) with 
X = 27 [28, 29] and 36 sites. Note that the X = 36 curve in Fig. H differs slightly 
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Figure 6. Magnetization curve of the 5 = 1/2 Heisenberg antiferromagnet on the 
Kagome lattice with N = 27 (dashed hne) and 36 sites (sohd line). The inset shows 
part of the Kagome lattice with a a/S x \/3 superstructure indicated by the circles in 
certain hexagons. 

from the corresponding curve in [29] for (M) < 1/3. The reason is that due to the 
computational effort only selected /c-sectors were investigated in [29] and in some cases 
the correct symmetry of the ground state was in fact missed. Our new = 36 curve 
agrees with unpublished results of Waldtmann and Everts [30] . 

The Kagome lattice is famous for a disordered ground state at h = accompanied 
by a small spin gap with many singlets inside this gap (see [31, 32] and references 
therein). The spin gap should give rise to a narrow (M) = plateau even if it is 
difficult to recognize in Fig. IHl 

A plateau at (M) = 1/3 can be recognized better in Fig. IHl In fact, the presence of 
this plateau at (M) = 1/3 in the 5 = 1/2 Heisenberg model on the Kagome lattice has 
been established previously by considering also system sizes different from those shown 
in Fig. ini [28,33]. Nevertheless, the state of this plateau is still an interesting issue. 
For the classical Heisenberg model at (M) = 1/3, thermal fluctuations select collinear 
states, but due to the huge degeneracy of this space, there appears to be no real order on 
the classical level at (M) = 1/3 [34] (see also [35]). For S = 1/2, it is useful to consider 
the XX Z model. In the Ising limit A ^ oo one can then first establish [36] a relation to 
a quantum dimer model on the hexagonal lattice which was argued [37, 38] to give rise 
to a valence bond crystal ground state with a \/3 x order indicated by the circles in 
the hexagons in the inset of Fig. IHl Furthermore, like for the triangular lattice one can 
compute the overlap (see Fig. El) of the (M) = 1/3 wave function of the XX Z model 
with that of the quantum dimer model as a function of A and one finds no evidence 
for a phase transition for A > 1 [36]. This implies that also the (M) = 1/3 state of 
the = 1/2 Heisenberg model on the Kagome lattice is a three-fold degenerate valence 
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bond crystal, i.e. an ordered state with many low-lying non-magnetic excitations, but 
still a small gap to a// excitations [36]. 

There may be a further plateau at (M) = 5/9 in Fig. IHl although it is difficult 
to draw unambiguous conclusions from the available numerical data in this region of 
magnetization values. 

Finally, one can see a pronounced jump of height 6{M) = 2/9 just below saturation 
in Fig.ini This jump arises due to local magnon excitations which can be constructed in 
the spin-S XX Z model on the Kagome lattice [29,39]. These local magnon excitations 
give rise to an exact degeneracy at the saturation field and thus a macroscopic jump in 
the magnetization curve of height 6{M) = 1/(95). In the state just below the jump, 
magnons are localized on the hexagons marked by circles in the inset of Fig. IHl Since 
this is a crystalline ordered state, general arguments (see e.g. [40,41]) predict a gap 
and consequently a plateau in the magnetization curve at (M) = 1 — 1/(95"). Indeed, 
a plateau at (M) = 7/9 can be seen clearly in the magnetization curve of the 5 = 1/2 
Heisenberg model on the Kagome lattice in particular when one considers also lattices 
with = 45 and 54 sites [29]. 

4. Conclusions and outlook 

In this article we have discussed examples of plateaus in one- and two-dimensional 
quantum spin models. In one dimension, the appearance of a plateau is controlled by 
the quantization condition Q which implies in particular that the magnetization (M) 
on a plateau must be rational. This condition can be interpreted as a commensurability 
condition for the ground state. Hence, it is important to observe that translational 
symmetry can be spontaneously broken in frustrated quantum spin models. Even in 
well-studied examples such as the frustrated 5=1/2 Heisenberg chain, an (M) = 1/3 
plateau with a spontaneously broken period three was discovered only recently [13]. 
Furthermore, an infinite sequence of plateaus in the orthogonal dimer chain [19] shows 
that there is no upper bound on the possible periodicity of spontaneous breaking of 
translational symmetry. 

It is less clear if the condition Q is also applicable to two and higher dimensions. 
Nevertheless, it works for most examples discussed here: The state of the (M) = 1/3 
plateau in the 5=1/2 triangular lattice has a unit cell with V = 3 spins (see section 
13. 2|) and for the Kagome lattice the plateaus with (M) = 1/3 in the 5=1/2 model 
as well as the one with (M) = 1 — 1/(95) in the spin-5 model both have unit cells 
with V = 9 spins (see section ESI in particular the inset of Fig. IHl). These cases fit well 
with the quantization condition 0. However, the disordered ground state at /i = in 
the 5 = 1/2 Kagome lattice (see [31,32] and references therein) is not covered by the 
condition Q since it gives rise to an (M) = plateau without the appropriate breaking 
of translational symmetry. 

From the point of view of experiments, the exchange constants Jij dictate the scale 
of magnetic fields needed to access the high-field region. SrCu2(B03)2 is an 5 = 1/2 
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system with sufficiently small J's that permit the observation of plateaus at (M) = 1/8, 
1/4 and 1/3 in pulsed-field magnetization experiments [42,43]. SrCu2 (603)2 is believed 
to be a good realization of the two-dimensional Shastry-Sutherland lattice [17]. Because 
of its relation to SrCu2(B03)2, this frustrated spin model has been analyzed theoretically 
in detail (see [44] for a recent review). We hope that the examples discussed in this 
article will stimulate further theoretical and experimental research on the magnetization 
process of frustrated quantum spin models. 
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